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Abstract. - We study the long-time asymptotical behavior of the survival probability P t of a 
tagged monomer of an infinitely long Rouse chain in presence of two fixed absorbing boundaries, 
placed at x = ±L. Mean-square displacement of a tagged monomer obeys X 2 (t) ~ t 1/2 at all 
times, which signifies that its dynamics is an anomalous diffusion process. Constructing lower and 
upper bounds on P t , which have the same time-dependence but slightly differ by numerical factors 
in the definition of the characteristic relaxation time, we show that Pt is a stretched-exponential 
function of time, In (Pt) ~ — t 1/,2 /L 2 . This implies that the distribution function of the first exit 
time from a fixed interval [— L, L] for such an anomalous diffusion has all moments. 



Introduction. — Dynamics of a tagged monomer 
(TM) of a polymer chain in solution is a practically im- 
portant physical example of an anomalous diffusive pro- 
cess. While dynamics of the whole chain is dominated 
by the motion of its center of mass which moves diffu- 
sively, dynamics of a TM is coupled to dynamics of other 
monomers; as time evolves, progressively more and more 
other monomers start to impede dynamics of the TM, 
which ultimately results in a subdiffusive motion. Such 
subdiffusive motion persists up to a certain characteris- 
tic time, proportional to some power of the chain length; 
at greater times conventional diffusive behavior is estab- 
lished [1]. 

For the so-called Rouse model of a polymer chain [2] , in 
which a chain is considered as a series of K beads linearly 
connected by harmonic springs, the TM mean-square dis- 
placement X 2 (t) can be calculated exactly. One finds an 
anomalous diffusion law X 2 (t) <~ \[t for times less than 
the so-called Rouse relaxation time Tr ~ K 2 , (the time 
needed for some perturbation, i.e. a kink, to spread diffu- 
sively along the whole chain), and conventional diffusive 
motion with reduced, by factor K, diffusion coefficient for 
times larger than Tr. In fact, for the Rouse model many 
important dynamical properties can be calculated exactly, 
e.g., the TM position probability distribution function, the 
dynamical structure factor [1], as well as the measure of 
different trajectories of a tagged monomer [3]. One may 
even determine exactly the dynamics of the TM of a Rouse 
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chain in more complex situations - in random layered 
flows [4, 5] or in situations appropriate to electrophoresis 
of polyampholytes, i.e., polymers whose monomers may be 
positively or negatively charged and the chain is subject 
to external electric field [6] . 

Recently, following a general interest in understanding 
subdiffusive and superdiffusive motion, a different aspect 
of anomalous diffusion of tagged monomers of a Rouse 
polymer chain has attracted some attention. Namely, dy- 
namics of a Rouse chain in presence of traps or absorbing 
boundaries reacting with some or just one of its monomers 
has been analysed [7-9]. Physically, such a situation is 
realized for polymers diffusing on solid surfaces contain- 
ing chemically active sites, which may react reversibly or 
irreversibly with any or some of the chain monomers tem- 
porarily or completely anchoring the chain. Conceptually, 
this question is interesting in its own right since the answer 
contains a solution, for anomalous diffusion, of a certain 
first-passage time problem, whose general understanding 
is a basic aspect of stochastic processes [10]. 

Dynamics of a TM of an infinitely long Rouse chain 
in one-dimensional systems in presence of two absorbing 
boundaries has been analysed in Ref. [8] within a path- 
integral formalism with an exact measure of trajectories 
of such a monomer [3], and a suitably extended classic 
method of images. It was shown that the probability Pt 
that the tagged monomer commencing its motion at the 
origin does not escape, during time t, from the interval 
[—L,L], or, in other words, that it "survives" up to time 
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t in presence of two traps placed at x 
obeys 

t l/2 



ln(P t )~- 



L 2 



-L and x = L, 



(1) 



i.e., the decay of Pt is described by a stretched-cxponcntial 
function of time. The law in Eq.(l) has been previously 
conjectured using heuristic arguments in Ref. [7]. 

On the other hand, dynamics of a tagged monomer of a 
finite Rouse chain in one-dimension in presence of absorb- 
ing boundaries has been analysed numerically in Ref. [9]. 
Here it was claimed that Pt is an exponential function of 
time 

ln(P t )~-i, (2) 

and thus decays at a faster rate. Consequently, the results 
of Refs. [7,8] and [9] are in an apparent contradiction with 
each other. 

In this Letter we aim to resolve this controversy by de- 
riving, using the approach outlined in Ref. [11], rigorous 
lower and upper bounds on the suvival probability Pt of 
a Rouse chain monomer in a one-dimensional system with 
two absorbing bondaries. We set out to show here that P t 
obeys the following double-sided inequality 
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In (P t ) /4L 



r 3/2 



< 1, 



(3) 



which defines the decay law up to a numerical factor in the 
characteristic relaxation time. This inequality confirms 
the result in Eq.(l) and rules out the result in Eq.(2). 

This paper is outlined as follows: In section 2 we present 
the notations and write down basic equations. In section 3 
we present the results of a heuristic approach, in which the 
survival probability of a tagged Rouse chain monomer in 
presence of two absorbing boundaries at x = ±L is inter- 
preted as the survival probability of a Brownian motion in 
presence of absorbing boundaries which move away from 
the origin as ztLt 1 / 4 . Next, in section 4 we derive rigorous 
lower and upper bounds on P t , which lead to the inequal- 
ity in Eq.(3). Finally, in section 5, we conclude with a 
brief recapitulation of our results and discussion. 

Notations and basic equations. Dynamics of 
a discrete Rouse chain comprising an infinite number of 
monomers is described by a set of Langevin equations [1] : 



dX n (t) 
dt 



= - (X n+1 (t) + X„_i(t) - 2X n (tj) + C t (n) , 



(4) 



X n (t) being an instantaneous position of the n-th bead, 
— oo < n < oo, and - independent Gaussian white- 
noise processes, such that 



d n) = o, 



^(n)^(m) 



Sn,mS(t - t'). 



(5) 



In Eq.(5), the bar denotes averaging over thermal histo- 
ries, S nim is the Kronecker symbol and 5(t) is the delta- 
function. Note that, for simplicity of presentation, we 
have set in Eq.(4) the friction constant equal to 1, the 



spring constant and the temperature equal to 1/2. These 
parameters can be easily restored in our final results. 

Note, as well, that Eq.(4) describes the time evolution 
of local heights of the Edwards- Wilkinson interface [12] 
in one dimension, or the time evolution of the difference 
of local concentrations of A and B species for diffusion- 
limited A + B — > inert reactions with random, steady, 
uncorrelated input of A and B [13]. Our results will thus 
apply to these systems too. 

Supposing that initially all monomers are at the origin, 
we have that X n= o(t) - position of the zeroth TM of an in- 
finitely long Rouse chain at time t - for a given realization 
of noises is determined as a portfolio of independent 
Gaussian processes: 

oo „ t 

X n=0 (t)=X(t)= ]T / dT(^e-^I n (t-T), 

n=-oo J 

(6) 

where I n (t — r) is the modified Bessel function of order n. 
Eq. (6) yields the following expression for the mean-square 
displacement of the zeroth monomer of an infinitely long 
Rouse chain: 



X 2 (t) =ter 2t [I (2t) + h(2t)] 



(7) 



As a matter of fact, any other initial condition can be 
considered. However, the effect of the initial state of the 
chain on the process in Eq.(6) fades out quite rapidly; it 
was observed in Ref. [9] that the difference in dynamics 
between the case when initially all monomers are at the 
origin or when one starts from an equilibrated configura- 
tion is rather small. Thus we have chosen the simplest 
case when X n (t = 0) = for any n. 

We note next that since we are concerned with the large- 
t behavior, it will not matter much how we define (| n ' - 
as continuous in time functions or as discrete processes, 
provided that we keep all essential features of noise. We 
thus divide, at fixed t, the interval [0,t] into N (N ^> 1) 
small subintervals A, (such that AN = t), and assume 



that Ct""* is constant and equal to (Y 11 /\/A within the 
fc-th subinterval, k — 0, 1, . . .,N — 1. We suppose that 
{Cfc"''} i s an infinite set of independent random variables 
with normal distribution N[0, 1]. 

Then, X(t) can be written down as a weighted sum of 
an infinite number of independent discrete noise processes: 



(n) 



oo N 



n= — oo 1=1 

with I and n-dependent weights 

r Al 



1 f 

<7* (n) = ^= / due"" !„(«). 
VA Ja(i-i) 



(8) 



(9) 



At this point, it is also expedient to introduce another 
property - an effective time-dependent variance of - which 
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will emerge in what follows as the key parameter. Squaring 
Eq. (8) and averaging the resulting expression with respect 
to distributions of i.i.d. variables {(^}, we get 



JV 



* 2 (*) = 5> 2 , 



(10) 



1=1 



for any k, the trajectory {Yk} is exactly the trajectory 
{Xk}, with the only difference that it evolves in the in- 
verse time N — k and is shifted by a constant (realization- 
dependent) value Xn. Consequently, the survival proba- 
bility P t = P(max\X(t)\ < L) = P(max|F(i)| < L). 

To calculate P(u\ax\Y{t)\ < L), we proceed as follows. 
Let I(max|y(£)| < L) be the following indicator function: 



where the effective variance of is given by 



oo „ 

*? = E fa B) ) = 

n— — oo 

- dm du 2 e- Ul - U2 1 (ui + u 2 ) . (11) 

A JA(l-l) JAIl-1) 



l(max\Y(t)\ <L) = 



1, max|F(i)| < L, 
0, max|F(£)| > L . 



(17) 



The integrations in Eq.(ll) can be performed exactly, but 
the resulting expression - a combination on nine modified 
Bessel functions - is rather cumbersome and is of a little 
use. In fact, all information we need to know about of 
can be extracted directly from Eq.(ll): 

(a) of is a monotonically decreasing function of I. To see 
this, it suffice to notice that exp{— x)Iq(x) is a monotoni- 
cally decreasing function of x; then, one finds from Eq.(ll) 
that of obeys the following double-sided inequality: 

Ae- 2Al I (2AZ) < of < Ae- 2A V-Vl (2A(Z - 1)) , (12) 

i.e., is bounded from both sides by monotonically decreas- 
ing functions of I. 

(b) d\ decays as when I — » oo. In fact, bounds in 
Eq.(12) become very sharp for I 3> 1 and 



In terms of descending partial sums Eq.(17) can be rewrit- 
ten as 

JV 

j(max|y(t)| <L) = l[l(\Y k \ <l). (18) 
fe=i 

Next, let recti (x) be a rectangular function, such that: 



recti (x) 



1, 

1/2, 
0, 



|x| < L, 
x = ±L, 
\x\ > L. 



(19) 



Representing recti(x) via its Fourier transform: 

dy sin(Ly) 



recti (x) 



J-oc K 



y 



exp [iyx] , 



(20) 



we write down the indicator function in Eq.(18) as the 
following TV-fold integral: 



~2 



1/2 



dy fe sin(Ly fe ) 



J/fc 



Inserting the latter expression into Eq.(10) and performing 
summation, we recover the result in Eq.(7). 

Define now the following event: An TV-step trajectory 
X(t), Eq.(8), commencing at the origin, docs not leave 
the interval [— L, L], L > 0, or, in other words, that 
max|JV(t)| < L. Such an event takes place, clearly, when 
the absolute value of any ascending partial sum 



x exp 



N 



E yfc Yk 



fc=l 



(21) 



Xk 



oo 



N 



/ j / j °l SAT-i! 

n= — oo l=N-k+l 



(14) 



Now, in order to determine the desired probability P t , 
we have to average the indicator function in Eq.(21) with 
respect to distributions of i.i.d. variables ■ To do this, 
we first rewrite the sum in the exponential in Eq.(21) in 
the following form 



JV 



JV 



which define positions of the tagged monomer at con- 
secutive discrete "time" moments k, k = 1,2, ... ,N, is 
bounded from above by L. 

However, in order to get a convenient "direction" of time 
in our final results, we will prefer to work with descending 
partial sums 



(22) 



n= — oo fc— 1 



m— k 



Inserting the latter expression into Eq.(21) and performing 
averaging, we find that P t is given by 



p, . r...rn* 

00 00 k=l 7 '" 



dy k sin(Ly fe ) 



Vk 



n— — 00 1=1 



"1 W-r 



Since, evidently, 



Y, 



N-k 



Xk = x 



N- 



(15) 



(16) 



x exp 



JV / JV 

E^ E 



fc=i 



i—k 



(23) 



where the effective variance of has been defined in Eq.(ll). 
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Next, changing the integration variables: 

Yi = Vi+V2 + ■■■+VN, 
Y 2 = V2 + V3 + ... + Vn, 
Y 3 = 2/3 + 2/4 + ... + VN, 



Y N = y N , 



(24) 



we obtain 



r°° f°° dY k 

* = J ■■■ II 

J — OO J —CO fc—i 



dY k sin (L(Y k -Y k+1 )) 



x exp 



--n 2 Y 2 



Y, 



N+l 



Y k - Y k+1 



0. 



(25) 



Now, it is expedient to use the following integral identity 
for the sine-function: 



sin(L (Y k -Y k+1 )) 
Y k - Y k+1 



i J dX k exp [iX k (r fe - Y k+1 )] . 

(26) 

Plugging Eq.(26) into Eq.(25), and performing integra- 
tions over Yfe-s, we finally arrive at the following meaning- 
ful representation of the survival probability: 



I : J :}\ 



dX k 



2ncr k 



exp 



(X k — X k _\) 



(27) 



where X n = 0. The integrand in Eq.(27) represents the 
measure of trajectories of the zeroth monomer of an infi- 
nite discrete Rouse chain. Its continuous-space counter- 
part has been calculated previously in Ref. [3]. 

Heuristic estimate of P t . — Changing the integra- 
tion variables, X k = <7 k x kl we can cast Eq.(27) into the 
following form: 



/ ■■■/ n 



V*s Jt dx . 



2tt 



x exp 
Further on, since 



1 / cr fc _i 

"7T I x k : — Xk-l 

Ok 



a k \k 



,x = 0. (28) 



(29) 



we may expect that, for sufficiently large k, this factor will 
not matter and 



x k 



-Xk-l 



(x k - Xk-i) ■ 



(30) 



Hence, Eq.(28) can be approximated by 

i-i/ffi r L/a 2 r L/a N N 



Pi 



-L/a\ J -L/5-2 J —L/o-n 
1 



dx k 



x exp 



(x k - x k -i 



2 



2ir 



x = 0. 



(31) 



One notices that the right-hand-side of Eq.(31) determines 
the probability that an TV-step Brownian motion trajec- 
tory does not leave an interval whose boundaries move 
dcterministically away from the origin as ±L/a k , i.e., P t 
can be approximately defined as 



/L/<j t 
P(X,t)dX, 
-LI a t 



(32) 



where a t - (A 2 /47rt) 1 / 4 , Eq.(13), while P{X,t) obeys 



dP(X,t) 1 d 2 P(X,t) 



dt 



2A dX 2 



P{X,t = 0) = S(X), (33) 



and 



P(X = ±L/at,t)=0. 



(34) 



We estimate next P t defined by Eq.(32) using an adiabatic 
approximation discussed in Ref. [15]. The basic idea be- 
hind this approximation is that, if the boundary advances 
sufficiently slowly, the density distribution approaches the 
same form as in the fixed boundary case, except that the 
parameters in this probability distribution acquire time 
dependence to satisfy moving boundary conditions [15]. 
We find that, within this approximation, in the leading 
order 



P(X, t) w exp 



( u: | ) (35) 



which yields the following estimate: 



exp 



7T 



1/2' 



(36) 



Note that this estimate agrees with Eq.(l). 



Rigorous bounds on the survival probability P t . 

— The TV- fold integral in Eq.(27) can not be, of course, 
performed exactly and recourse has to be made to control- 
lable approximations. Below we construct rigorous lower 
and upper bounds on P t in Eq.(27), which both have the 
same time dependence defining in such a way an asymptot- 
ically exact (up to a constant factors in the characteristic 
relaxation time) result. 

The method we use here is based on the approach out- 
lined in Ref. [11] within the context of the Riemann- 
Liouvillc fractional Brownian motion. In constructing 
bounds, we will take advantage of the following two facts: 

(a) The effective dispersion a k in Eq.(27) is a monotoni- 
cally decreasing function of time. 

(b) A fundamental property of P t defined by Eq.(27) is 
that P t — Pti&i, <?2, CT3, . . . , ctjv) is a monotonically de- 
creasing function of any variable 5 k [11]. 

This signifies that replacing any or all a k by £(fc), such 
that <7fc < S(fc), we will decrease the survival probability 
and arrive at the lower bound on P t ; if, on contrary, we 
will replace one or all a k by S(fc), such that a k > S(fc), we 
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will increase the survival probability and obtain an upper 
bound on P t . 

We start with an upper bound on P t . Since u k is a 
monotonically decreasing function of k, we have that 

a k > t(k) = a N , (37) 

which holds for any k. The equality is attained only for 
k = N. 

Hence, Pt in Eq.(27) is bounded from above by 



We notice now that the expression in the second line in 
Eq.(45) defines the probability that an TV-step trajectory 
of Brownian motion, evolving in time T t , Eq.(44), does 
not leave the interval [-L(JL.)V4 ; L(^) 1 / 4 ]. Hence, P t 
in Eq.(27) is bounded from below by 



Pt > exp 



7T 

8L 2 



4A 



1/2 



— T, 



exp 



7T 



1/2 



L N 



Pt < 



/Li pLi 
/ 'I 



dX k 



L/ctjv 

-L/5 N 



k=l 
L/a N 



exp 



(X k — X k _i) 



N 

n 

0. 



dxi 



exp 



(x k - x fc _i) z 



where Xo = and xq = 

Expression in the second line in Eq.(38) determines the 
probability that an TV-step Brownian motion trajectory 
does not leave an interval with fixed boundaries ±L/<jn, 
which is a classic problem in the probability theory (see, 
e.g., Ref. [14]). Consequently, we find that at sufficiently 
large times 

1/2 

= exp 



(46) 

This bound rules out an exponentially fast decay of the 
survival probability suggested in Ref. [9]. 

Finally, combining the lower and the upper bounds, we 
(38)>btain the double-sided inequality obeyed by P t , Eq.(3), 
which defines the decay of the survival probability P t up 
to a numerical factor in the characteristic relaxation time. 



Pt < exp 



- — 

8 L 2 



16L 2 



(39) 



Consider now a lower bound on P t in Eq.(27). To con- 
struct such a bound, we turn back to the definition of 
<r k in Eq.(ll). Recollecting that exp(— x)Io(x) is a mono- 
tonically decreasing function of x, we have that for any 
u\ 1 u 2 > 0, 

e- Ul - U2 I (u 1 +u 2 )<e- U2 I Q (u 2 ). 
Consequently, b\ is bounded as 

5%< alu 2 e- U2 h (u 2 ) . 

JA(fc-l) 

Further on, since for any u 2 > 0, 

e- U2 Io(u 2 ) < -±=, 



(40) 



(41) 



we have that the following bound holds: 



~o\ < E 2 (fc) = 2 



A 



1/2 



(T fe — T fe _i) , 



(42) 



(43) 



where 

T k = k 1 ' 2 . (44) 

Hence, the survival probability P t in Eq.(27) is bounded 
from below by 



L N 



> 



/Li f'Li ly 

I II 



dX,, 



H- 



-.1/4 



k=l 



2?rS(fc) 

1/4 n 



exp 



(X k — Xk-i) 
2E 2 (fc) 



x exp 



(Xh - Xh-l) 



2 



n 



dxi 



2 (T k — 7fe_i) 



x = 0,x = 



Conclusions. — To conclude, we have studied the 
long-time asymptotical behavior of the probability Pt that 
a tagged monomer of an infinitely long Rouse chain will 
not escape, up to time t, from an interval [—L,L]. We 
have shown, by constructing rigorous lower and upper 
bounds on P t , which both have the same dependence on 
time but slightly differ, by numerical factors, in the defini- 
tion of the characteristic relaxation time, that Pt follows 
ln(P t ) ~ — t 1 / 2 /^ 2 . This decay law confirms our earlier 
predictions based on uncontrollable approaches [7, 8] and 
contradicts a recent prediction ln(P t ) ~ —t based on nu- 
merical simulations [9]. 

This result implies that the probability distribution 
function of the first-exit time from an interval [-L, L] for 
the anomalous diffusion process executed by the tagged 
monomer of an infinitely long Rouse chain has a stretched- 
exponential tail <~ exp(— t x l 2 jL 2 ) and thus has all mo- 
ments. 

We note that the obtained decay law agrees, as well, 
with the general result on the survival probability of the 
Riemann-Liouville fractional Brownian motion (fBm) in 
presence of absorbing boundaries [11] when the Hurst in- 
dex H = 1/4. To realize that this is not a coincidence 
and the dynamics of a tagged monomer of an infinitely 
long Rouse chain is indeed in the fBm " universality class" , 
consider Eq.(4), in which, for simplicity, we treat n as a 
continuous variable. Then, we have that X n =o(t) obeys 



0. 



(45) 



X n=0 (t) 



Crdr 



(t-r) 



1/2 



dn exp 



2(t-r) 



In 



(47) 

where £t is a white noise in time and /„ is a white 
noise of variable n, — oo < n < oo. Notice now that 
exp[— n 2 /2(t — t)] is a bell-shaped function which broad- 
ens with time, which signifies that as time progresses, 
more and more monomers start to affect the dynamics 
of the zeroth monomer. Let us replace, for simplicity, 
exp[— n 2 /2(t -— t)] by a rectangular function rect(n), such 
that rect(n) = 1 for n < \Jt — t and rect(n) = for 
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n > y/t — t. Then, Eq.(47) reads [14] Feller W., An Introduction to Probability Theory and 

Its Applications (Wiley, New York) 1971 
/"* C T dT f^ 1 ^ [15] Krapivsky P.L. and Redner. S., Am. J. Phys., 64 (1996) 

*n=o(*)« / ^72 / _dnf n . (48) 546 . 

JO {t — Tj ' J-Vt^T 

Now, one expects that, for typical realizations of /„, 

[ A dnf n ~ A 1 / 2 , (49) 
and hence, X n=0 (t) follows 

x n=0 (t)« f-^-m, (50) 

Jo (t - t) 7 

which is exactly the Riemann-Liouville fractional Brown- 
ian motion with Hurst index = 1/4 studied in Ref. [11]. 
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